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SUMMARY 


The  high  cost  of  modern  weapons  systems  forces  increasingly 
stringent  requirements  on  structural  integrity  and  durability 
of  new  military  hardware.  Prediction  of  deterioration  of  material 
properties  in  storage  and  service  environments  requires  therefore 
an  increase  in  sophistication  which  necessitates  a strong  interaction 
between  experiment  and  analysis. 

This  report  is  Part  IV  of  an  investigation  into  the 
moisture  effects  in  carbon  fiber  reinforced  epoxy  composites. 

Specif ically, it  deals  with  the  prediction  of  changes  in  the  elastic 
behavior  of  these  composites  as  a result  of  moisture  uptake  in  a 
natural  environment. 

This  progr^  was  funded  by  the  Naval  Air  Systems  Command 
Task  No.  A320000'0010123  during  the  period  of  May  1,  to 
December  3,  1976^' 

J.  R.  DIXON 
By  direction 
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INTRODUCTION 

This  report  is  part  IV  of  an  investigation  into  the  effects 
that  moisture  has  on  carbon  fiber- reinforced  epoxy  composites. 

Part  I discussed  the  diffusion  of  moisture  in  composites^  part  II 
summarized  our  experimental  observations  on  the  reversible  changes  in 
strengths  of  composites  after  moisture  exposure,  part  III  described 
an  analytical  modeling  of  uptake  and  internal  distribution  of  mois- 
ture in  composites  with  respect  to  real  outdoor  weathering  conditions. 
This  report  discusses  the  changes  in  the  elastic  behavior  of  composite 
as  a result  of  outdoor  exposure. 

A unified  scheme  for  the  prediction  of  the  elastic  behavior 
of  multi-angleplied  laminates  is  pretented.  It  includes  temperature 
effects,  swelling  effects  and  changes  in  the  state  of  stress  (strain) 
of  a laminated  composite  as  a function  of  time  under  any  prescribed 
climatic  environment. 


DISCUSSION 


A.  Problem 


It  has  been  observed  by  various  investigators  that 

exposure  of  fiber-reinforced  epoxy  composites  to  moisture  leads  to 
a reduction  in  matrix  dominated  strength  properties.  The  degree  of 

CT3  J.  Hertz,  Final  Report  NASA  8-27435  Jmrs  1973  (NASA  Contract) 
C23  c.  E.  Browning,  28th  Annual  Conference  of  the  Society  of 

Plastic  industry,  Feb  1973  (Washington,  D.C.)  Proceedings  15A. 
C3l  J.  M.  Augl  and  A.  E.  Berger,  "Moisture  Effects  on  Carbon  Fiber 
Epoxy  Composites".  Proceedings  p.  SAMPE  Meeting,  Seattle, 
Washington,  Oct  1976. 

C4J  J.  M.  Augl,  "The  Effect  of  Moisture  on  Carbon  Fiber-Reinforced 
Epoxy  Composites.  II  Mechanical  Property  Changes." 

NSWC/WOL/TR  76-149  (1976) . 

C5l  D.  H.  Kaelble,  P.  J.  Dynes  and  L.  Maus,  "Hydrothermal  Aging  of 
Composite  Materials  Part  1:  Interfacial  Aspects."  Journal  of 
Adhesion  8121  (1976)  Part  2:  Matrix  Aspects  ibid  8,  155  (1976) 
C63  eTlT  McKague,  J.  E.  Halkias  and  J.  D.  Reynolds  "Moisture  in 
Composites.  The  Effect  of  Supersonic  Service  on  Diffusion." 

J.  Comp.  Mat.  9,  2 (1975). 

Vf]  Proceedings  Air  Force  Workshop  on  Durability  Characteristics 
of  Resin  Matrix  Composities.  Battelle  Columbus  Laboratories, 
Sept  30,  1975. 
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strength  reduction  depends,  of  course,  on  the  type  of  failure 
mechanism,  on  the  moisture  concentration,  and  on  the  temperature. 

The  strength  degradation  is  reversible,  i.e.,  after  removal  of  the 
absorbed  moisture  the  original  strength  can  be  restored.  Also,  the 
strength  loss  levels  out  after  the  equilibrium  concentration  with 
the  environment  has  been  reached  £43. 

In  Part  I of  this  report  £8)  we  have  discussed  the  diffusion  of 
moisture  in  composites  and  in  Part  III  £93  we  have  described  an 
analysis  of  how  real  outdoor  and  service  conditions  affect  the  uptake 
and  the  internal  moisture  distribution  in  composites. 

In  the  following  we  shall  discuss  how  a real  outdoor  and  service 
environment  affects  the  elastic  properties  of  carbon  fiber  reinforced 
epoxy  composites,  i.e.,  the  change  of  the  elastic  properties  as  a 
function  of  time  in  a prescribed  climatic  environment. 

We  shall  restrict  ourself  to  laminated  plates  that  are  composed 
of  transversely  isotropic  unidirectional  plies  layed  up  at  different 
fiber-angles.  We  also  assume  that  the  plate  is  exposed  to  moisture 
on  both  surfaces. 


An  example  of  such  a composite  is  shown  in  Figure  1 which  gives 
a 4 ply,  +450,  symmetric  laminate.  The  curved  line  (shaped  like  a 
parabola)  indicates  a superposed  profile  of  moisture  concentration 
after  a certain  time  of  moisture  exposure  which  is  given  by  the 
solution  of  the  diffusion  equation  as  discussed  in  reference  183  and 
D3.  The  moisture  distribution  is  symmetric,  however,  the  individual 
plies  are  no  longer  transversely  isotropic  nor  are  they  symmetric 
because  the  elastic  moduli  of  the  matrix  is  not  constant  anymore  in  the 
2-direction  since  the  modulus  of  the  matrix  is  dependent  on  the 
moisture  concentration  £43.  We  can,  therefore,  no  longer  take  the 
individual  elastic  constants  for  the  laminae  that  constitute  the 
laminate  but  we  have  to  know  the  change  of  the  lamina  elastic  constants 
as  a function  of  absorbed  moisture,  and,  if  we  are  interested  in  the 
temperature  effects,  we  have  also  to  determine  the  lamina  elastic 
constants  as  a function  of  temperature. 


However,  we  can  approximate  this  continuous  change  in  the  matrix 
modulus  in  the  z-direction  by  finite,  discontinuous,  steps  (A) 
"sufficiently  small"  to  neglect  the  concentration  differences  on  either 
side  of  the  A-layers,  and  we  consider  these  A-layers  again  as  having 
transversely  isotropic  symmetry.  What  is  "sufficiently  small"  is  a 
question  of  the  change  in  modulus  with  moisture  concentration  and  the 
sensitivity  of  such  measurements.  The  A differences  should  depend  on 
the  gradient  of  the  moisture  concentration,  and  are  equal  or  smaller 
than  the  physical  plies. 


£83 

£83 


J.  M.  Augl  and  A.  E.  Berger,  "The  Effect  of  Moisture  on  Carbon 
Fiber- Reinforced  Composites.  I Diffusion."  NSWC/WOL/TR  76-7 
(1976)  . 

J.  M.  Augl  and  A.  E.  Berger,  "The  Effect  of  Moisture  on  Carbon 

Fiber-Reinforced  Composites;  III  Prediction  of  Moisture 

Sorption  in  a Real  Outdoor  Environment."  NSWC/WOL/TR  77-13  (1977), 
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Thus,  the  whole  problem  reduces  again  to  classical  laminated 
plate  theory,  which  has  been  amply  described  in  various  papers  and 
monographs  of  which  we  want  to  cite  only  a few  |l0  - 14J. 

This  theory  may  be  expressed  by  the  compact  matrix  equation 

(1) 

l]-[lj§][|-]  '!> 

where  the  N matrix  is  the  stress  resultant  matrix;  M is  the  moment 
resultant  matrix;  e°  is  the  matrix  describing  the  midplane  stretching 
and  shearing  deformations,  and  < describes  twisting  and  bending 
deflections;  the  A,  B,  and  D matrices  are  the  extensional  stiffness 
matrix,  the  bending-stretching  coupling  matrix  and  the  flexural 
stiffness  matrix  respectively. 

The  following  steps  are  therefore  required: 

1.  Determine  the  moisture  diffusion  coefficient  of  the  composite 
in  the  z-direction  as  a function  of  temperature  (and,  if  necessary, 

as  a function  of  concentration) . 

2.  Calculate  the  internal  moisture  concentration  distribution 
for  a prescribed  environment  (with  variable  temperature  and  humidity) 
as  a function  of  time  by  solving  Pick's  second  diffusion  equation. 

3.  Determine  the  constitutive  properties  Qij  of  each  A - layer 
as  a function  of  uniform  moisture  distribution  and  as  a function 

of  temperature. 

4.  Calculate  the  stiffness  matrix  of  the  laminate  (A,  B,  D)  from 
the  corresponding  A - 

We  have  discussed  step  1 and  2 previously  bf  9]  . In  the 
following  we  shall  investigate  step  3 and  step  4. 


tPJ-  J.  E.  Ashton,  J.  C.  Halpin  and  P.  H.  Petit,  "Primer  on  Composite 
Materials."  Composite  Materials  Workshop,  Technomic  Publ. 

_ Co.,  Inc.,  Stanford,  Conn.  (1969) . 

J.  C.  Halpin  and  S.  W.  Tsai,  "Environmental  Factors  in  Composite 
Materials  Design"  AFML  TR  67-423. 
jl^  J.  R.  Vinson  and  T.  W.  Chou,  "Composite  Materials  and  their 
Use  in  Structures"  John  Wiley  1975. 

[13)  S.  W.  Tsai,  "Mechanics  of  Composite  Materials.  Part  II." 

AFML-TR  149  (1966). 

^ly  S.  G.  Lekhnitskii,  "Anisotropic  Plates."  Translated  by 

S.  W.  Tsai  and  T.  Cheron,  Gordon  and  Breach,  1968. 
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B.  A-Layer  Constitutive  Properties 

For  an  orthotropic  lamina  oriented  in  the  principal  material 
direction,  such  that  the  fibers  are  aligned  parallel  to  the  1-axis  as 
shown  by  the  sketch  below 


23 

'13 

12 


'll 

^12 

^13 


constitutive 

relations 

are 

Qi2 

Qi3 

0 

0 
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- 

^1 

^22 

Q23 

0 

0 

0 

^2 

Q23 

®33 
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2Q44 
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**^23 

0 

0 

0 

2Q55 

0 

0 

0 

0 

0 

2066 

_;*Y12^ 

where  O2,  o^,  "^22'  ^13'  ^12  normal  and  shear  stresses, 

^1'  ^2'  ^3'  ^23'  "'^13'  ''^12  corresponding  strains  and  the  . 

are  the  components  of  the  stiffness  matrix.  ^ 
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We  have  used  the  momenclature  of  reference  ClO,  12l . Since  we 
can  assume  that  the  fibers  are  parallel  to  the  1 axis  but  randomly 
distributed  in  the  plane  perpendicular  to  the  1 axis,  this  2-3-plane 
is  a plane  of  isotropy  and  the  lamina  is  said  to  be  "transversely 
isotropic"  or  sometimes  also  called  "monotropic."  For  such  a system 
the  following  relations  hold: 


^33  ®22i®13  °12-.°55  ®66 

and  = Js  (Q22  - Q23) 


Further,  in  thin  plate  and  shell  theory,  we  assume  that  there  are  no 
stresses  applied  in  the  z-direction  (=  3-direction)  therefore 

®13  " ®23  ^ ®33  " 

A detailed  discussion  of  the  elastic  stiffness  matrices  for  all 
the  crystallographic  symmetry  classes  was  given  by  S.  G.  Leckhnitskii 
ri53  (Attention  should  be  given  to  the  confusion  that  may  arise  from 
naming  the  axes  1,  2,  3,  and  z,  y,  z by  different  authors  with  respect 
to  fiber  direction  and  plane  of  isotropy.) 

The  physical  meaning  of  the  matrix  elements  are  as  follows: 


o = 

^11 

°11 

(1-V12  V22I 

0 = 

®22 

^12 

(I-V22  V21> 

Qi2  = 

'’21  _ "12  ®22 

21  ^21^  ^^“^12  ^21^ 

Q44 

S3 

Q55  = 

S6  " ^12 

Where  En 

and  E22  the  Young's  moduli  in 

[15]  S.  C.  Lekhnitskii,  Theory  of  Elasticity  of  Anisotropic 
Elastic  Bodies.  Holden  Day  1963. 
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If  the  stiffnesses  of  a plate  are  to  be  determined  for  some  new 
reference  axes  (x,  y)  in  the  plane  of  the  lamina,  this  can  be 
accomplished  by  a matrix  transformation  describing  a rotation  of  the 
coordinates  around  the  3-axis  (=  z-axis)  by  an  angle*  0 as  indicated 
by  6-8. 


H(xy) 

H (XY) 


and  thus 


[o] 


(xy) 


where 


It]. 


2 2,, 

m n 0 

2 2 „ 

n m 0 


0 0 
0 0 
-mn  mn 
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[t]-i  . 


m 

n2 

0 

0 

mn 


and  where  m s cos 


= H(12| 


' It]-1  [Q.  .]  [t] 


0 2mn 
0 -2mn 
m-n  0 
n m 0 
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~2mn 
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2mn 
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0 

0 

-n 
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0 

mn 

0 

0 

/ 2 2, 
(m  -n  ) 

9 , n = sin  9 . 


(6) 

(7) 

(8) 


{■ 

i 


i 


*The  angle  is  (+)  if  it  is  measured  counterclockwise  from  the 
axis  to  the  (+1)  axis. 


w 


NSWC/WOL  TR  77-61 

Thus  a stiffness  matrix  (Qij)  is  obtained  describing  a generally 
orthotropic  lamina.  Stress  strain  relations  for  such  a generally 
orthotropic  lamina  including  transverse  shear  deformation  are  given 
by  (9) 


where  [q  ] = [t]'!  [q  ] [t] 


yz  (9) 


The  evaluation  of  the  Qij  was  simplified  by  Tsai  and  Pagano  ri6i 
who  used  certain  functions  that  are  independent  of  0 and  are  given  by 
the  following  relations 


®22  ” ^1  ~ ^2  ^ ^3 

Qi2  = - U3  cos  (46) 

Qgg  = cos  (40) 

= -^5  U2  sin  (20)  - U3  sin  (40) 
Q2g  = -*5  U2  sin  (20)  + U3  sin  (40) 


where 


“1  = 8 Ou  * 3 O22  + 20^2  * ^066> 

'’2  ' ■ ®22’ 

= I^Qll  ^ Q22  - 2Qi2  - 

^4  ^ ^ ^22  ^ ^^12  " ^^66^ 

U = i 

5 ®(Qii  O22  ■ ^®12 

S.  W.  Tsai  and  N.  Pagano,  "Invariant  Properties  of  Composite 
Materials."  Composite  Materials  Workshop.  Technomic  Publ. 
Co.,  Inc.,  pp.  213-53  (196S). 
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and  the  quantities  Q.  . Qj.^  Q*c  were  given  by  J.  R.  Vinson  and 
T.  W.  Chou  [123, 


Q44 

0 

It 

2 

cos 

9 

2 

+ Qgg  sin 

9 

Q55 

0 

II 

. 2 
sin 

9 

2 

+ Qgg  cos 

9 

(12) 

Q45  = {Q55  - Q44)  cos  9 sin  9 


The  determination  of  the  required  A-layer  constitutive  properties 
can  be  either  measured  directly  on  unidirectional  composite  specimens 
or  predicted  from  the  constituent  properties. 

a.  Direct  measurement  of  A-layer  constitutive  properties. 

To  obtain  all  the  elastic  constants  of  a lamina  for  various 
levels  of  moisture  and  temperatures  is  a rather  involved  experimental 
program. 

In  order  to  prepare  the  test  specimens  it  is  necessary  to  expose 
the  samples  at  various  constant  relative  humidity  levels  for  a long 
enough  time  to  reach  a uniform  moisture  concentration.  For  a .25  cm 
thick  composite  this  may  take  aboui.  80  days  even  at  750c.  The 
following  measurements  of  composite  elastic  constants  will  have  to  be 
carried  out  at  about  5 different  temperatures  (between  room 
temperature  and  170oc)  and  5 different  humidity  levels:  Eiif  ^22'  ^12’ 
V12f  G23*  If  only  a minimum  of  5 specimens  are  used  per  condition, 
this  would  amount  to  625  specimens. 

b.  Prediction,  of  A-Layer  Constitutive  Properties  from  the  resin 
and  fiber  properties. 

Numerical  techniques  have  been  developed  to  predict  the  composite 
elastic  constants  from  the  known  properties  of  the  matrix  and  fiber. 
These  problems  can  be  solved  rigorously  if  a periodic,  geometric, 
array  of  fibers  is  assumed  (for  instance,  tetragonal  or  hexagonal). 
Also,  for  a statistical  distribution  of  fibers,  solutions  have  been 
described . 

Even  if  the  prediction  of  properties  is  not  always  very  accurate, 
these  formulas  are  quite  useful  for  our  purpose  of  predicting  changes 
of  properties. 

The  reason  why  the  theories  do  not  always  predict  accurately  the 
experimentally  determined  values  may  be  due  to  the  fact  that  one  or 
several  of  the  assumptions  made  for  the  model  are  not  correct.  Most  of 
these  theories  make  the  following  assumptions: 

1.  The  ply  is  macroscopically  homogeneous,  linearly  elastic,  and 
transversely  isotropic. 
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2.  The  fibers  are  linearly  elastic  and  homogeneous. 

3.  The  matrix  is  linearly  elastic  and  homogeneous. 

4.  Fiber  and  matrix  are  free  of  voids. 

5.  There  is  complete  bonding  at  the  interface  between  fiber  and 
matrix  and  there  is  no  transition  region. 

6.  The  fiber  and  matrix  are  initially  in  a stress  free  state. 

7.  The  fibers  are  regularly  spaced  and  aligned. 

In  a real  composite  neither  assumption  is  totally  justified. 
Composites  have  voids,  the  fibers  are  not  uniform,  they  are  not 
regularly  spaced,  they  are  not  perfectly  aligned,  there  is  not 
everywhere  a perfect  bound  between  resin  and  fiber,  there  is 
experimental  (microscopic)  evidence  that  there  is  a transition  zone 
between  resin  and  fiber,  simple  thermal  calculation  show  that  there  may 
be  considerable  residual  stresses  in  composites  (with  high  Tg  resins) 
that  have  been  cured  at  a high  temperature.  It  is  therefore  rather 
surprising  how  well  the  models  predict  the  composite  properties. 

What  we  are  interested  in  here  is  not  so  much  the  absolute  values 
of  the  elastic  constants  but  rather  the  changes  in  properties. 

Even,  if  the  assumptions  are  true  only  to  a certain  extent, the 
violation  will  remain  relatively  constant  if  we  let  moisture  permeate 
into  the  resin.  There  will  be  no  change  in  the  number  of  voids;  the 
relative  fiber  spacings  remains  constant;  the  misalignment  remains 
unchanged;  the  transient  zone  remains,  although  it  may  change 
properties;  the  interfacial  bonding  may  be  weakened,  this  is 
particularly  true  for  glass  fiber  composites;  moisture,  however,  may 
enhance  relaxation  of  residual  stresses. 

Thus,  after  the  initial  composite  unidirectional  properties  have 
been  experimentally  determined  in  the  dry  state,  as  baseline 
properties,  we  expect  that  the  relative  changes  in  the  elastic 
composite  properties,  due  to  moisture  sorption,  can  be  predicted  with 
a considerably  higher  degree  of  accuracy  than  it  can  be  done  from 
original  resin  and  fiber  constituent  properties,  so  that  it  is  only 
necessary  to  determine  the  matrix  elastic  constants  Em,  Gm  (or  vm)  as 
a function  of  uniform  moisture  concentration  and  temperature. 

These  properties  can  be  determined  nondestructively  on  very  few 
simple  plate  specimens  by  ultrasonic  sound  velocity  measurements  of 
longitudinal  and  shear  waves.  Torsional  pendulum  experiments  combined 
with  a few  tensile  measurements  may  also  be  used.  Since  there  is  a 
certain  moisture  desorption  at  the  high  temperature  measurements,  it  is 
preferable  to  use  not  too  thin  plate  specimens  (.3  cm  thickness  will 
probably  suffice) . 
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We  made  no  attempt  here  to  review  the  various  approaches  that  have 
been  taken  to  predict  the  elastic  constants  of  unidirectional 
composites  but  rather  list  in  Appendix  A a number  of  solutions  for  the 
different  elastic  constants  that  have  proved  to  be  quite  useful  in 
predicting  these  properties. 

C.  Properties  of  a Laminate  with  Absorbed  Moisture. 


After  moisture  has  penetrated  the  composite,  the  modulus  of  the 
matrix  is  not  only  changed  in  the  neighborhood  of  the  glass  transition 
temperature  (Tg)  but  over  a very  wide  temperature  range  (see  Figure  2) . 

Since  the  matrix  modulus  enters  in  all  the  prediction  formulas  of 
the  composite  elastic  constants  (see  Appendix  A) , it  is  obvious  that 
the  composite  properties  change  accordingly.  If  we  know  from 
experiment  the  degree  of  change  of  the  matrix  with  moisture  and 
temperature  and  if  we  know  the  internal  moisture  distribution,  we  can 
describe  in  a straight  forward  fashion  the  composite  elastic  behavior, 
by  using  the  same  formalism  as  is  used  for  analyzing  laminates  with  a 
homogeneous  matrix  and  different  ply-orientations.  Instead  of  treating 
the  laminate  as  being  composed  of  real  plies  we  describe  it  in  terms  of 
the  above  mentioned  artificial  A-layers,  and  we  simply  rewrite  the 
equations  of  laminated  plate  theory  in  terms  of  A-layers  as  if  the 
composite  were  composed  of  a hybrids  of  A-layers.  In  the  following  we 
adapt  the  basic  nomenclature  used  in  reference  flO,  12}. 


Hki  =[o]k4  I'l 


where  a and  e are  the  stress  and  strain  matices  and  Q is  the 
stiffness  matrix  with  respect  to  the  laminate  axes  (x,y,z),  the 
subscript  kA  refers  to  the  kth  A-layer.  This  equation  can  be 
rewritten  in  terms  of  strain-displacement  relations 


LOkA  = ^ 


where 


e°  = — 

X 3x  i 
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substituting  (14)  into  (13)  gives 


WkA  = [^IkA  ^ ^ [q]  kA  [^1 


The  stress  and  momentum  resultants  can  be  written  in  the  usual  way  by 


h/2  a 


T I dz,  I M 

xy  y 


-h/2  T, 


h/2  a 


-h/2  a 


a I zdz 

y 


where  the  N N N are  the  in-plane  stress  resultants,  M^,  M „ 
X y xy  2 

are  the  stress  couples.  Qx  and  Qy  are  the  transverse  shear  forces 

and  h is  the  thickness  of  the  laminate. 

Summation  over  all  A-layers  yields 


V y h ' 

^ i ^A  -1  ["xy 


C / WkA  dz  + y- 

= ^ ^A-1  e°  . ^A-1 


X I zdz 


lince  and  [<lare  not  functions  of  z,  one  can  write  (17)  as 


kA  kA 

S l"°l  f <5"  + [O]  k,  (kl  / 


‘kA-1 
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or  in  short  form 

[n]  = [a]  [e°]  + [b]  [k]  (19) 

where 

n 

A.  • = (Q.  (h,,.  - h,  . (20) 

13  / ^ '^1]  kA  kA  kA-1)  ' ' 

kA  = 1 

n 

«i3  ' E (21) 

kA  = 1 

(i,j  = 1,2,6) 


Similarly  for 
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if  the  transverse  shear  forces  are  to  be  included  the  assumption 
is  made  that  the  transverse  shear  stresses  are  distributed  paraboli- 
cally  across  the  laminate  thickness  and  a weighing  function  as  used 
by  Reissner  employed 


the  transverse  shear  resultants  are  obtained  for  the  k A - layer 


(kA)  ^ Q(kA)  1 ^2  iy 

^xz  ^55  2 ^xz  ^^45  2^yz 


^ ^45  2'xz  ^44  2^yz 


and  from  (16)  one  obtains 


where 


^ 2^  xz  ^ ^45  2^yz^ 

""  ^^^45  2^  xz  ^44  2^yz^ 


- 4 


4 ^^kA  " 


1 3 


(i,j  = 4,5) 


Finally  the  inclusion  of  all  artificial  A-layers  leaves  the  overall 
formalism  unchanged  and  it  can  be  written  as  is  customary  by 


All  ^12  ^^16  I ®11  ®12  ^®16 


^12  ^22  ^^26  j ®12  ®22  ^®26 

^16  ^26  ^^66  I ®16  ®26  ^®66 


1 

IIU. 


®11  ®12  ^®16  j °11  °12  ^°16 

®12  ®22  ^®26  I °12  °22  ^'°26 

®16  ®26  ^®66  j °16  °26  ^°66 


y 

1 

2^xy 


E.  Reissner,  "On  a Variational  Theorem  in  Elasticity.' 
J.  Math,  and  Phys.  29,  90  (1950). 
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and  by  (26),  or  by  the  abbreviated  formula  (1)  which  comprise  the 
full  system  of  constitutive  equations  to  study  laminates  that  are 
symmetrically  permeated  by  moisture.  Other  useful  relations  between 
N,  M,  A,  B and  D are  given  in  Appendix  B. 

The  changes  with  time  are  governed  by  the  diffusion  equation. 

D.  Property  Changes  in  Laminates  due  to  the  Effect  of  Temperature 
On  the  Matrix  Modulus 

The  matrix  modulus  changes  with  temperature  and,  increasingly  so 
as  the  Tg  of  the  matrix  is  approached. 

Again,  the  formalism  discussed  in  the  previous  sections  remains 
the  same. 

The  Pick's  diffusion  equation  has  the  same  form  also  for  thermal 
diffusion.  The  only  difference  is  that  the  thermal  diffusivity  is 
several  orders  of  magnitude  larger  than  the  moisture  diffusion.  The 
thermal  diffusion  is  so  much  faster  that  for  most  problems  in 
composite  structures  we  are  more  interested  in  the  isothermal  steady 
state  effects  at  different  temperatures,  or  with  constant  temperature 
gradient  rather  than  the  transient  phenomena  of  non-steady  state  heat 
flux.  Though  there  is  no  difficulty  in  dealing  with  transient 
conditions  because  of  the  difference  in  rates  of  thermal  and  moisture 
diffusion.  The  changes  in  the  elastic  behavior  due  to  temperature 
changes  and  transien*-  internal  heat  fluxes  is  simply  superposed  on 
the  "constant"  internal  moisture  distributions.  Even  the  same 
A-layer  spacings  may  be  used.  The  changes  of  the  properties  with 
time  are  now  governed  by  the  thermal  diffusion  equation,  which 
requires  that  the  thermal  diffusivity  of  the  composite  in  the 
z-direction  be  known. 

The  matrix  elements  of  the  stiffness  tensor  (Qij)kA  replaced 
by  the  elements  (Qij  (CH20rT))kA  indicating  that  the  A layer  elastic 
constants  depend  on  internal  moisture  concentration  and  temperature. 

E.  Property  Changes  in  a Laminate  due  to  Volume  changes  caused 
by  solvent  swelling  or  temperature  changes. 

In  addition  to  affecting  the  modulus  of  the  matrix,  moisture 
sorption  and  temperature  changes  have  similar  effects  in  that  they 
change  the  specific  volumina  of  the  constituent  materials  unequally. 
This  leads  to  a change  in  the  internal  state  of  stress. 

We  do  not  intend  to  discuss  this  subject  here  further  since  it 
was  well  described  by  Pipes,  Vinson  and  Chou  in  a recent  paper  C187. 


R.  B.  Pipes,  J.  R.  Vinson  and  T.  W.  Chou,  "On  the  Hygrothermal 
Response  of  Laminated  Composite  Systems."  J . Comp . 

Materials  10,  129  (1976). 
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These  authors  have  modified  the  constitutive  equations  of  the 
laminated  plate  theory  to  include  these  effects.  And,  for  a given 
laminate  strain,  the  stresses  within  each  lamina  were  given  by 


Oj  '^(z,t)  = Q^jJ^ej°(t)  + zKj  (t) -aj'^T(z,t) M(z,t)J 
(i,j  = 1,2,6) 


(29) 


where  a j , Qjj#  and  icj  are  the  stresses,  stiffness  components, 

inplane  strains  and  bending  strains  respectively;  aj  and  6j  are  the 
thermal  and  hygroscopic  expansion  coefficients  and  T(z,t)  and  M(z,t) 
are  the  temperature  and  moisture  concentration  distributions,  at  the 
time  t.  The  superscript  k refers  to  the  ply. 


This  formula  still  neglects  the  temperature  and  humidity 
dependence  of  the  matrix  modulus. 


F.  Combined  Effects  of  Temperature  and  Moisture  on  the  Elastic 


Behavior  of  Composites 


Expanding  on  the  effort  of  Pipes,  Vinson  and  Chou  we  propose  to 
modify  their  equation  (29)  to  include  the  changes  in  the  matrix 
modulus  and  therefore  in  the  A-layer  stiffness  components  to 
the  following  form: 


I”’'’'’  So/'Jks  ^ “"j"’  - 


,(‘'3  • 1.2.6) 


(30) 


where  the  symbols  are  the  same  as  in  equation  29  , however  the 
(Qij)'s  are  now  functions  of  temperature,  moisture  concentration  and 
time,  whereby  the  concentration  of  moisture,  CH20f  and  the  time  t 
are  related  by  the  Pick's  diffusion  equation,  and  the  temperature 
and  time  are  similarly  related  by  the  analogous  equation  for  heat 
conduction.  However,  because  of  the  large  rate  differences  in 
thermal  and  moisture  diffusion,  it  is  more  convenient  to  treat  these 
two  phenomena  separately. 


G. 


Computer  Program  for  Calculating  the  Elastic  Stiffness  Matrix 
of  a Laminate  as  a Function  of  a Real  Outdoor  Environment 


For  calculating  the  change  in  the  elastic  behavior  of  a composite 
for  real  outdoor  environment  we  require  the  following  input  data: 
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1.  Long-term  average  temperature  and  relative  humidity  of  the 
environment,  or  1-3  hourly  weather  data  of  past  years  on  magnetic 
tape. 

2.  Diffusion  coefficient  of  the  laminate  in  the  z-direction 

as  a function  of  temperature  (and  concentration,  if  necessary) , or  the 
diffusion  coefficient  of  the  matrix  as  a function  of  temperature  (and 
concentration)  together  with  the  composite  fiber  volume  fraction. 

3.  Equilibrium  concentration  of  the  matrix  material  as  a 
function  of  relative  humidity  (and  temperature,  if  necessary). 

4.  Fiber  and  matrix  moduli  (Em  and  Gm  or  Vj^,  Enf,  E22f»  Gl2f» 
vi2f)  as  a function  of  uniform  moisture  concentration  and  temperatue. 
(The  fiber  properties  are  considered  to  remain  constant  over  the 
temperature  range  of  interest.)  For  predicting  the  lamina  properties 
Ell  Ei2f  Gi2  and  vi2  the  Halpin-Tsai  equations  (A17-A19)  were  used  and 
for  G23  Rosen's  Equation  (A16)  was  used. 

5.  Laminate  geometry  (number  of  plies,  orientation  of  plies, 
thickness  of  plies,  thickness  of  A-layers) . 

6.  Time  intervals  of  printout  of  the  stiffness  matrix. 

7.  Applied  laminate  stresses  or  strains  (aj,ej). 

The  printout  for  the  specified  time  intervals  gives: 

1.  The  (if  desired) 

2.  The  matrices  A,  B,  D,  A*,  B*,  C*,  D* 

k k 

3.  The  strains  and  stresses  in  each  lamina  (e.;  , o-  ) or,  if 

desired  of  each  A-layer.  ^ ^ 

H.  Examples. 

Example  1:  Composite  stiffness  changes  as  a function  of  matrix 
stiffness  changes. 

Moisture  will  plasticize  the  matrix  and  reduce  the  values  of  the 
matrix  elastic  properties  depending  on  the  amount  of  (uniformly 
distributed)  moisture  and  on  the  temperature. 

Let  us  assume  that  these  values  have  been  determined 
experimentally,  what  is  the  corresponding  change  in  the  lamina 
elastic  constants. 

Table  1 gives  the  original  fiber  matrix  and  lamina  properties. 
Table  2-6  and  Figures  3-7  show  the  calculated  changes  in  the 
unidirectional  lamina  properties  as  a function  of  change  in  the 
matrix  modulus  for  various  fiber  volume  fractions  (V^) . 
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Table  1.  Elastic  Moduli  of  Narmco  5208/T300  graphite  Fiber  lamina 
and  its  constituents.  (The  values  are  given  in  10^  psi.) 


Elast.  Const, 


Narmco  5208 
Epoxy 


T300  Composite 

Carbon  Fiber  5208/T300  .62V^ 


^11 

.5 

33.5 

20.8 

^22 

.5 

3.05 

1.46 

^12 

.185 

3.0 

.633 

^12 

.35 

.285 

.31 

Example  2: 

Quick  Estimate 
Composite. 

of  the 

Ultimate  Stiffness 

Change 

We  would  like  to  know  how  long  it  takes  for  a 4-ply 
unidirectional  T300/5208  composite  to  reach  equilibrium  in  an 
environment  such  as  the  Panama  Canal  Zone,  and  what  would  be  the 
expected  percentage  reduction  of  the  elasticconstants  of  such  a 
composite  when  tested  at  100°C  (212°F) , as  compared  to  the  dry 
composite  properties  tested  at  the  same  temperature. 

Known  data; 

Average  temperature  in  the  Canal  Zone  is  about  77°F  and  the 
humidity  if  fairly  constant  about  93%  RH. 

The  thickness  of  a ply  of  Narmco  5208/T300  after  fabrication 
into  a laminate  is  about  .01397  cm. 

The  diffusion  coefficient  of  the  5208/T300  composite  at  this 
temperature  is  about  1.0  x 10“^0  cnr/sec. 

This  problem  is  very  simple  and  can  be  solved  without  further 
use  of  a computer. 

From  reference  8,  Figure  4,  we  know  that  a 95%  saturation  value 
is  given  by  Dt/h2  = 1.25  (where  D is  the  diffusion  coefficient, 

h is  ^ the  thickness  of  the  laminate  and  t is  the  time  to  reach 

this  percentage  of  equilibrium) . 

Since  the  thickness  of  a 4-ply  laminate  is  .05588  cm,  the  time 
to  reach  95%  of  the  equilibrium  concentration  is  about  16  weeks  (see 
Figure  17  of  reference  8). 
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Table  2.  Calculated  change  of  the  lamina  Young's  modulus,  En, 
of  a T300  composite  as  the  matrix  modulus  Em  changes.  (Values 
are  given  in  106  psi.) 


.4 

.5 

.6 

.7 

.5 

13.70 

17.00 

20.30 

23.60 

.4 

13.64 

16.95 

20.26 

23.57 

i 

.3 

13.58 

16.90 

20.22 

23.54 

1 

13.52 

16.85 

20.18 

23.51 

.1 

13.46 

16.80 

20.14 

23.48 

Table  3.  Calculated  change  of  the  lamina  transverse  Young's 


modulus,  E22r  of  a T 
(Values  are  given  in 


.5 

.4 

.3 

.2 

.1 


00  composite  as  the 
106  psi) 


.4 

.5 

1.005 

1.189 

.856 

1.029 

.688 

.844 

.496 

.622 

.271 

.349 

matrix  modulus  changes. 


. 6 

.7 

1.411 

1.682 

1.244 

1.515 

1.042 

1.304 

.790 

1.023 

.459 

.623 
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3C 


Table  4.  Calculated  change  of  the  lamina  longitudinal  shear 
modulus,  Gi2,  of  a T300  composite  as  the  matrix  shear  modulus, 
Gm/  changes.  (Values  are  given  in  106  psi.) 


; \ f 
m \ 

.4 

.5 

.6 

.7 

200 

.415 

.511 

.642 

.832 

175 

.368 

.455 

.576 

.753 

150 

.320 

.399 

.506 

.668 

125 

.271 

.338 

.433 

.577 

100 

.220 

.276 

.356 

.479 

075 

.167 

.211 

.274 

.374 

050 

.113 

.144 

.188 

.260 

025 

.057 

.073 

.097 

.135 

Table  5.  Calculated  change  of  the  lamina  transverse  shear 
modulus,  G23/  of  a T300  composite  as  the  matrix  shear  modulus 
Gm  changes.  (Values  are  given  in  106  psi.) 


iVf 

m \ 

.4 

.5 

.6 

.7 

.2 

.882 

.879 

.946 

1.085 

.175 

.809 

.804 

.867 

1.000 

.150 

.730 

.723 

.781 

.907 

.125 

.641 

.634 

.687 

.803 

.100 

.543 

.535 

.581 

.685 

.075 

.433 

.425 

.463 

.552 

.050 

.308 

.302 

.330 

.397 

.025 

.166 

.161 

.177 

.216 

» 
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Table  6.  Calculated  change  in  the  lamina  major  Poisson's  ratio, 
vi2f  of  a T300  composite  as  the  matrix  Poisson's  ratio  Vni  changes. 
(Values  are  given  in  106  psi.) 


i vVf 

! m ' 

.4 

.5 

.6 

.7 

.5 

.414 

.392 

.371 

.349 

.4 

.354 

.342 

.331 

.319 

.3 

.294 

.292 

.291 

.289 

.2 

.234 

.242 

.251 

.259 

.1 

.174 

.192 

.211 

.229 
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From  Figure  2 we  see  that  at  100°c  the  retention  in  the  matrix 
modulus  is  81%  of  the  room  temperature  value  for  the  dry  resin  and 
only  41.5%  for  a resin  equilibrated  at  93%  RH,  i.e.,  if  the  room 
temperature  resin  modulus  was  .5  x 10^  psi>  the  100°C  moduli  for  dry 
and  wet  resin  are  .405  x 10°  and  .207  x 10°  psi. 

From  Figures  3-7  we  can  read  off  directly  the  changes  in  elastic 
constants  of  the  composite.  Assuming  a fiber  volume  fraction  of  .6 
we  see  from  Figure  3 that  no  measurable  change  in  En  is  expected; 
from  Figure  4 one  can  see  that  E22  will  drop  from  1.46  x 10°  to 
1.25  X 10^  and  .8  x 10®  psi  respectively;  from  Figure  5 one  can  see 
the  longitudinal  shear  modulus  G12  will  drop  from  .633  x 10^  psi 
to  .512  X 10°  and  .263  x 10°  psi  respectively;  from  Figure  6 we  see 
that  the  transverse  shear  modulus  G23  will  drop  from  .90  x 10°  to 
.729  X 10°  and  .373  x 10®  psi  respectively. 

A compar ison^of  the  moduli  of  the  dry  and  moisture  exposed 
composites  at  100  C shows  that  there  is  a 50%  retention  in  the  matrix 
modulus,  no  change  in  the  tensile  modulus  En,  a 64%  retention  in 
the  transverse  tensile  modulus,  a 51%  retention  of  the  longitudinal 
modulus  Gi2  and  a 51%  retention  of  the  transverse  shear  modulus  G23. 
Since  it  is  not  expected  that  the  Poisson's  ratio  of  the  resin  will 
change  significantly  below  glass  transition  temperature,  the  change 
in  the  composite  Poisson's  ratio  will  be  even  less  (see  Figure  7). 

The  implications  are,  that  not  only  are  certain  stiffness 
parameters  reduced,  but  also  the  failure  strength  will  be  reduced 
if  the  strain  to  failure  is  not  increasing  correspondingly  in  order 
to  permit  the  same  state  of  internal  stress. 

At  the  present  time  we  do  not  have  enough  information  about 
the  change  in  strain  to  failure  of  the  matrix  material  as  a function 
of  temperature  and  absorbed  moisture  to  compare  experimental  and 
calculated  values,  therefore,  such  considerations  are  still 
qualitative. 

Example  3;  The  Temporal  Change  in  Elastic  Properties  of  a composite 
under  Washington  D.  C.  Weather  Conditions. 

An  8 ply  unidirectional  and  an  8 ply  +45°  symmetric  laminate 
(5208/T300;  Vf  = .7)  are  thought  to  be  exposed  undei,  Washington  D.  C. 
conditions.  The  properties  of  the  matrix,  fibers  and  dry  composte 
are  given  in  Table  1. 

The  internal  distributions  after  0.1,2,3,5,7,8,28,45  and  50 
months  were  calculated  as  indicated  in  reference  9.  From  Figure  2 
(and  Table  1 and  7)  and  by  applying  the  A-layer  analysis  the  A,  B 
and  D matrices  were  calculated  for  the  composite  to  be  tested  at 
100°C.  Only  the  A matrices  for  the  exposure  times  of  0,3,8  and 
45  months  exposure  are  given  in  Table  8.  Figure  8 shows  the  change 
in  Age  (the  shear  stiffness  of  the  8 ply  unidirectional  laminate). 


i 

i 
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Table  8.  "A  Matrices"  for  8 -ply  T300/Narmco  5208  (Vj  - . 
laminates  (0  - 0*  and  0 - ±45*  symm.)  after  various 
exposure  times  to  Washington,  D.C.  weather  conditions 
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It  also  shows  the  change  in  strain  the  8 ply  +45°  symmetric  composite 
would  show  under  a load  of  500  pounds/inch  in  the  0°  direction  during 
this  exposure  (creep  deformation  excluded).  One  can  see  that  these 
changes  are  proportional  to  the  amount  of  moisture  uptake. 

Example  4:  Effect  of  Internal  Moisture  Distribution. 

It  is  of  interest  to  answer  the  question  of  how  important  it 
is  to  know  the  internal  moisture  distribution  in  a composite  in 
addition  to  the  moisture  content  (the  latter  could  be  determined 
nondestructively) . If  there  is  an  unsymmetric  moisture  distribution 
in  a symmetric  composite  then  the  B matrix  is  no  longer  zero  but 
if  the  moisture  distribution  is  symmetric  the  B matrix  remains  zero. 

We  have  chosen  6 cases; 


Case  1:  The  same  composite  as  given  in  the  previous  example. 


Case  2 to  Case  6 are  shown  in  Figure  9 which  show  composites 
with  arbitrary  symmetric  ply  orientation  and  symmetric 
distribution  of  moisture,  with  unsymmetric  ply  orien- 
tation but  symmetric  moisture  distribution,  and  with 
unsymmetric  ply  and  moisture  distributions.  For  the 
cases  2-6  a non  linear  Em  vs  H2O  concentration  was 
assumed  of  the  form  + e2  = l.  The  Poisson's  ratio 
of  the  matrix  was  assumed  to  remain  constant  (v  = .35), 
the  fiber  proper  ties  are  those  of  T300  and 

°n,  = ^ 

Table  8 shows  the  difference  in  the  composite  of  case  1 with  the 
actual  and  with  a uniform  moisture  distribution  of  the  same  overall 
concentration.  In  this  case  there  is  no  significant  difference. 

The  A matrices  of  the  cases  2-6  are  given  in  Table  9;  also  given 
are  the  percent  differences  in  strains  between  the  composites  with 
uniform  and  nonuniform  moisture  distributions  if  a load  of  500  lbs 
would  be  applied  in  the  zero  direction. 

It  is  obvious  that  in  general  there  is  a difference  in  the 
elastic  behavior  between  composites  with  nonuniform  and  uniform 
moisture  distribution  but  with  the  same  overall  moisture  content. 

It  is  not  surprising  that  it  shows  up  most  strongly  in  simple 
+450  composites  (case  3 and  case  6)  though,  even  there,  with  a 
^natural"  distribution  obtained  under  real  conditions  as  in  case  1, 
the  % difference  in  strain  (at  a load  of  500  Ibs/inch)  does  not 
exceed  1%  during  the  entire  exposure  time.  Thus,  in  a first 
approximation  one  can  say  that  the  effect  of  the  internal  moisture 
distribution  obtained  under  real  conditions  may  be  neglected 
compared  to  the  stronger  effect  of  the  overall  moisture  content. 


Table  9.  Effect  of  Internal  Moisture  Distribution 
on  the  Elastic  Behavior  of  Composites. 
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CONCLUSIONS 

A compact  formalism  has  been  proposed  to  predict  the  changes 
in  the  elastic  behavior  of  fiber  reinforced  composites  by  combining 
real  outdoor  weather  data  and  a small  number  of  laboratory  measure- 
ments. The  analysis  is  based  on  a combination  of  the  theoretical 
concepts  of  diffusion,  micromechanics  and  a modified  laminated  plate 
theory.  It  was  also  demonstrated  that,  while  the  overall  moisture 
content  in  composites  can  change  significantly  the  elastic  properties, 
the  changes  of  elastic  properties  due  to  the  internal  distribution 
obtained  under  natural  conditions  does  not  differ  much  from  that 
of  a uniform  distribution  with  the  same  overall  moisture  content. 

More  work  is  required  for  predicting  the  change  in  composite 
strength  as  a result  of  moisture  penetration. 

ACKNOWLEDGEMENT 

The  authors  would  like  to  thank  Messrs.  M.  Stander  and 
C.  Bersch  of  the  Naval  Air  Systems  Command  for  their  interest  and 
support  of  this  work. 


25 


DNOD  O^-H 


NSWC/WOL  TR  77-61 


FIG  1 SYMMETRIC  MOISTURE  DISTRIBUTIONS  IN  A FOUR-PLY  COMPOSITE. 


E^.  RESIN  MODULUS  (10®  PSD 


FIG. 


CHANGE  OF  COMPOSITE  MODULUS  (E22)  AS  A FUNCTION  OF 
RESIN  MODULUS  IN  T 300  CF  UNIDIRECTIONAL  COMPOSITES 


COMPOSITE  SHEAR  MODULUS  (G12)  (lO^  PSI) 
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0.2  1.5  1.0  0.5  0 

Gm,  RESIN  MODULUS  (lOOpsi) 

FIG.  5 CHANGE  OF  COMPOSITE  SHEAR  MODULUS  (G12)  AS  A FUNCTION  OF 
RESIN  SHEAR  MODULUS  IN  T300CF  UNIDIRECTIONAL  COMPOSITES 
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CHANGE  IN  SHEAR  STIFFNESS,  STRAIN  AT  CONSTANT  LOAD  (500  LB). 
AND  MOISTURE  UPTAKE  OF  AN  8 PLY  LAMINATE  5208/T300  AS  A 
FUNCTION  OF  TIME  (WASHINGTON  DC) 
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APPENDIX  A 

PREDICTIONS  OF  THE  ELASTIC  CONSTANTS  OF  UNIDIRECTIONAL 
COMPOSITES  FROM  THE  CONSTITUENT  PROPERTIES 

Longitudinal  Modulus  E^^^  (=  Young's  Modulus) 

The  longitudinal  modulus  is  usually  described  by  the  rule  of 
mixtures 


E 


11 


V, 


+ E V 
m m 


(Al) 


or  modified  to  include  a fiber  misalignment  factor  k 


Ell  " ^^®f  '^f  ^ 


m m 


(A2) 


the  rule  of  mixture  was  further  modified 


by  Ekval  iri9j 


E 


11 


E,  V, 


+ E'  V 
m m 


E' 


E 

m 


l-2v 


2 

m 


where  Ef  is  the  longitudinal  fiber  modulus,  Em  is  the  matrix 
modulus,  Vf  and  Vm  are  the  fiber  and  resin  volume  fractions  and 
Vm  is  the  resin  Poisson's  ratio. 


(A3) 


Transverse  Modulus  E22 


Haskin  and  Rosen  derived  the  following  relation  ^20, 


Li9l  J.  C.  Ekvall,  AIAA  6th  Structural  and  Materials  Conf.  Palm  Springs, 
_ Apr  1965. 

L20J  Z.  Haskin  and  B.  W.  Rosen,  "The  Elastic  Moduli  of  Fiber- 
Reinforced  Materials."  J.  Appl.  Mech.  Jun  1964,  p.  223. 

L21J  B.  W.  Rosen,  N.  F.  Don  and  Z.  Haskin  "Mechanical  Properties 
of  Fiber  Composites."  NASA  CR-31  (1964). 


A-1 
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22  ^^33  K 


^^23^23 


23 


iliG 


23 


4)  = 1 + 


^^23  ^12^ 
^11 


K23  and  G23  are  bulk  and  shear  modulus  governing  plane-strain 
deformation  in  the  2,3-plane. 


Tsai  proposed  [223 


K,(2K  +G  )-G  (K.-K  )V 


m m 


W ^ /V 

f m m 


+ C 


K.  (2K  +G  .)  +G,  (K  -K,) 
t m t t m f m 

(2K  +G,)-2(K  -K,)V„ 
m t m r m 


] 


where 


Kj  = Ej/2(l-v^) 

K = E /2(l-v  ) 
m m m 

G^  = E^/2(1+Vj) 

G = E /2(l+v  ) 
m m m 


* fiber  Poisson's  ratio 

= matrix  Poisson's  ratio 
m 

c = o for  isolated  fibers 
c = 1 for  fibers  in  contact 


or  approximated  by 


E^  E 
f m 


'22 


E V.  + E^V 
m f f m 


Major  Poisson's  ratio  Vj^2 


= V-V  + V V 
12  f f mm 


(=  rule  of  mixture) 


|[223  S.  W.  Tsai,  "Structural  Behavior  of  Composite  Materials," 
NASA  CR-71  (1964)  . 


(A4) 


(A5) 


(A6) 


(A7) 
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A refinement  by  Rosen  £20^ 


E,V.L,  + V E L»v 
- f f 1 m m 2 n> 

12  ~ V,E,L,  + V E L- 
f f 3 mm2 


(A8) 


L = 2v  (l-v;)V  + V^(l-v„)v 

-i  £ m £ m m : 

= V^(l-v^-2v^) 

L = 2(l-v„^)V.  + (l+v„)V 
j m £ mm 


m 


Tsai  propc‘?ed  (A9)  [22^ 

K-V,(2K  -G  )V,+K  V (2K.+G  )V 
- ff  mmfmm  fm  m 

'"12  ■ ^ rrK~G  1 ^ tk^-k  ■) V 

f m m m f m m 


+ C 


+G.)V  +K  v.(2K  +G,)V. 
mm  £ £ m t £ m £ £ 

K,  (2k  +G.)  +G,  (K  -K,)  V_ 


t in  r f in  f in 
where  the  symbols  are  the  same  as  defined  under  (A5) . 
Longitudinal  Shear  Modulus,  Gj^2 


(A9) 


Ekvall  ^23^  proposed  (AlO)  and  (All) 


G,G 
f m 


'12  V G,  + V,G„ 
m £ £ m 


(for  square  filaments)  (AlO) 


where 


G G 

’12  " rFg"  r'(T-R)G;  filaments) 

m £ 


(All) 


2Tt 


-f 


sin  Ode 


(G  /G,)  + sin  0 1-(G  /G,) 
m £ m £ 


4V, 


R = 


4V,  4(1-V^) 

£ _ m 


V (V,  + V ) TT  TT 

£ m 


[23]  7.  C.  Ekvall,  ASME  paper  61  AV-56,  Aviation  Conf.,  Los  Angles, 

Mar  1961. 


A- 3 
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Rosen  t20']  proposed  an  upper  (AlO)  and  a lower  bound  (All) 


for  G 


12 


G,  _ (upper  bound)  =»  G (m_V, +V_) 

-L  z in  G i ^ 

G,-  (lower  bound)  = G (V,/m_)  + V_ 
Jt  ^ m X G ^ 


(A12) 

(A13) 


where 


V.  = 

1 2/3 


Vj  = 1-V^ 


rn 


: n(H-3^)  + (1-n^) 

n(l-(5^)  + (1-3^) 


yy  - 

n ^ 


» ' v, 


for  a random  array  he  proposed 


n(l+V,)  + 

c*  s X ^ m 

^12  m Tiv“+r+v:~~~ 
m f 


Tsdi  £22]  proposed 


(A14) 


2Gf-(G.-G^)V  (G,+G_) -(G,-G  )V 

p — M — r*  X z _ rn  , t m x m m 

^12  m 2G  +(G.-GJV"  ^ TGr+GT+rGl-GTv" 

m fmm  fm  fm  m 


(AIS) 


Transverse  Shear  Modulus  G 


23 


Rosen  proposed  £20,  2l] 


[(a-V^)  (l+pV^)-3Vj(l-V^)^6^  3 


(A16) 


wher :? 


a - 


P 


Y-1'' 


^ ” *\  iT 


m 


1 

3-4  V 


!ni_” 

l+yS 


y6^ 


^f  ~ 3-4'v, 
m f 


'^F 


= G„  “ Gp 

n2  ^23 


A- 4 
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Halpin-Tsai  Equations. 


For  design  purposes  and  rapid  computational  procedures  Halpin 
and  Tsai  11  have  proposed  a set  of  equations  which  are  simpler 
than  most  of  the  above  formulas  but  are  good  enough  for  estimating 
with  reasonable  accuracy  the  ply  properties; 


wi\er  e 


, = E,V,  + E V 

11  f f mm 

V,  - = + V V 

12  f f mm 


? (l+CnVj) 

'n,  ' fT-n-v->- 


[(Pf/P_)-l]/  [(p  /p 


p = composite  moduli  ^22/  ^12'  ^23' 

= corresponding  fiber  modulus  *^12f  ^12f  ^’23f^  ' 

p " corresponding  matrix  modulus  E , G , v ; 
ir.  ^ ^ m m m' 

r,  ■-  3 measure  of  reinforcement  which  depends  on  the  boundary 
conditions.  Approximate  values  are  given  by  = 2(a/b)  and 

^22 

by  log  rG, , - 3 log  (a/b)  , i.e.,  for  round  fibers  Cr.  - 2 and 

'’12 


{A17) 

(A18) 


(A19) 


APPENDIX  B 

Useful  relationships  between  A,  B and  D matices. 

Prom 
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APPENDIX  C 

PROGRAM  TO  CALCULATE  LAMINATE  MATERIAL  PROPERTIES 

Below  we  outline  a computer  program  which  calculates  the  time 
history  of  the  elastic  properties  of  a laminate  comprised  of 
transversely  isotropic  unidirectional  plys  of  a carbon  fiber 
reinforced  epoxy  composite  when  both  faces  of  the  laminate  are 
exposed  to  prescribed  climatic  conditions. 

A program  is  already  available  which  calculates  the  time 
dependent  moisture  distribution  within  such  a laminate  C3  , 9 1 . The 
computer  program  described  in  this  section  calculates  the  elastic 
properties  in  such  a laminate  when  given;  its  geometry,  its  internal 
moisture  distribution,  its  temperature,  the  dependence  of  the  elastic 
properties  of  the  fiber  and  of  the  resin  on  internal  moisture 
concentration  and  temperature,  and  the  dependence  of  the  elastic 
properties  of  a ply  of  the  composite  upon  the  elastic  properties  of 
the  fiber  and  of  the  resin. 


NOTATION 

The  laminate  lies  in  the  region  -z  < z ^ z.  This  region  is 
sub'l:  /ided  into  A-layers  whose  endpoints  = P(k)  are  located  at 

-z  = P(l)  < P(2)  < ...  < P{NDL+1)  = z 


where  NDL  stands  for  number  of  delta  layers.  These  points  will  be 
symmetrically  located  about  0 = P(1  + NDL/2) , and  each  interval 
[^P(k)f  P(k  + i)^  will  lie  within  one  ply  of  the  laminate. 

The  values  of  internal  moisture  concentration  (at  any  given  time 
t)  at  the  points  P(k)  determined  using  the  program  described  in  ([93 
are  denoted  by  C(k)  (k=l,  ...,  NDL+1) . The  points  P(k)  are  chosen 
to  be  spaced  closely  enough  so  that,  as  far  as  determining  the 
material  properties  within  each  delta  layer  Ak  = CP(k),  P(k+l‘/3  the 
moisture  concentration  within  Ai<  can  be  taken  to  have  the  constant 
value  W = (C(k)  + C(k+l)/2.  Assuming  that  the  elastic  properties  of 
each  delta  layer  can  be  calculated  given  W and  the  temperature,  the 
problem  of  determining  the  elastic  properties  of  the  laminate  is 
reduced  to  classical  laminated  plate  theory,  as  described  earlier  in 
this  text. 
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FORMULAS  USED  FOR  LAMINA  ELASTIC  CONSTANTS 

The  matrix  elastic  constants  Em  and  (or  Gja  = 2^/(2  (l+Vm) ) 
are  assumed  to  be  given  by  formulas  or  tables  as  functions  of  uniform 
internal  moisture  concentration  and  temperature,  and  the  fiber 
elastic  constants  Ehf,  Gi2Ff  vi2f  are  taken  to  remain  constant. 

For  predicting  the  lamina  properties  En,  £22^  Gl2»  and  vi2»  t^e 
Halpin-Tsai  equations  (A17-A19)  were  used  (here  the  1-direction  is 
parallel  to  the  fibers) ; 


®11  = ^llf  Vf  ^ ^m  ^m 


'’12  “ ''l2f  '^f  '^m  '^m 


(Cl) 


(C2) 


where 

fraction, 


= fiber  volume  fraction  and  V^,  = 1-V,:  = matrix  volume 


Also 


(14cvV,) 

P = ^ p 

(l-nV.)  *^m 


(C3) 


where 


V = 


[(Pf/P^)-!]  / t(P,/P„)  t 0 


and  where  P is  either  E22  or  G]^2  with  Pf  the  corresponding  fiber 
modulus  (E22f  or  Gi2f)  and  with  P^  the  corresponding  matrix  modulus 
(Eip  or  Gm) . For  round  fibers  the  parameter  ^ is  given  by  ?(E22) 

= 2 and  ^ (G12)  = 1.  For  calculating  G23,  Rosen's  equation  (A16) 
was  used; 


®23  = ^m 


(a+Pm^f)  (l+PVf)-3V^(l-V^)8^ 
(a-Vf) (1+pvJ)-3Vj(1-Vj)2b2 


m 


(C4) 


where 


Y+3 


' * 3-1 V ' “ 


m 


m 


3-4Vj'  P 


1+Y6f' 


Y * 


*12  f 


m 


C-2 
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We  now  give  the  notation  used  in  the  computer  program  for  the 
quantities  in  equations  (10),  (12),  (15),  (16),  (27),  and  (28). 

The  angle  of  the  k^h  delta  layer  (the  counterclockwise  angle  in 
degrees  through  which  the  positive  x-axis  must  be  rotated  to  reach 
the  positive  1-axis)  is  denoted  by  8(k).  For  each  k;  Qn,  Qi2/  022/ 
044/  055/  066  are  given  by  (3)-(5); 


'11 


'11 


<^“'^12''21 


) ' ®22 


'22 


(I-V12V21)'  ®44  “ ®23' 


Q55  “ °66  ®12' 


'12 


'21 


^21  ^11 
<l-^12'’2l' 


^12  ^22 
(I-V12V21) 


Def  i 

Lne  m = cos(9(k))  and  n = sin(6(k)) 

and  define 

oil 

BAR(k)  = 

Oil 

m + 2(Qi2+2Qgg) 

m2n2 

m n 

■f 

Q22 

BAR(k)  » 

OiiH 

+ 2(Qi2+2Qgg)m 

2 2 

+ 

^22” 

Q12 

BAR(k)  = 

<Qii 

'*’^22~^®66^ m + 

4 4 

O12 (m  +n  ) 

Q66 

BAR(k)  = 

(Oil 

+Q22-2Qi2-2Q66)"‘ 

2n2  + 

4 4 

®66  ^ 

Q16 

BAR(k)  = 

(Qii 

(O12- 

Q22‘'2Q66^'‘'^'" 

Q26 

BAR(k)  = 

(Oil 

-0i2~2Q66) n m + 

(O12- 

®22''’^®66^ 

Q44 

BAR(k)  = 

Q44 

m^  + Q55n^ 

Q55 

BAR(k)  = 

Q44 

n^  + Q55m^ 

Q45 

BAR(k)  = 

(055 

-Q44) mn 

For 

ifj  = 1/ 

2,6  define: 

NDL 

''ij 

= r 

BAR(k) • (P(k+1)- 

P(k)) 

k=l 

NDL 


B.  . = 

13 


^ Q. j BAR(k)* (P(k+l)‘-P(k)‘)/2 


(C5) 


(C7) 


k=l 


NDL 

Dij  « Q..  BAR(k)*  (P(k+l)^-P(k)^)/3 


(C8) 


C-3 


The  3 by  3 matrix  A is  defined  to  have  the  entries 

(All  A12  2A16\ 

A12  A22  2A26  ) 

A16  A26  2A66/ 

and  similarly  with  B and  D. 

The  3 by  3 matrices  ASTAR,  BSTAR,  CSTAR,  DSTAR  are  given  (in 
terms  of  matrix  algebra)  by  (Appendix  B) 

ASTAR  = a“^ 

BSTAR  = -A~^B  (C9) 

CSTAR  = BA“^ 

DSTAR  = D-BA"^B. 
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